## Alternant matrix

from PyM import *

[_,a,b,c,d,e] = multivariate_polynomial_ring(Z_,'a','b','c','d','e')

V = vandermonde([a,b,c,d,e],4)

show(V)

p = 19; K = Zn(p)

a = [element(j,K) for j in range(2,p) if is_prime(j)]

V_1 = vandermonde(a,5,1)

h = [element(j,K) for j in range(1,1+len(a))]

A = alternant_matrix(h,a,5)

show(A)